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ABSTRACT 

Electric field autocorrelations for the two-dimensional electro- 
static guiding-center plasma are calculated numerically. It is con- 
cluded that the autocorrelation, averaged over a thermal equilibrium 
ensemble, is damped in an approximately exponential fashion, as pre- 
dicted by Taylor and McNamara. Oscillatory behavior of the type 
predicted by Taylor and Thompson is not observed. 
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I. INTRODUCTION 


A central quantity in the theory of the transverse diffusion 

in a two-dimensional electrostatic guiding-center plasma 1 ” 10 is the 

-¥ 

autocorrelation of the kth Fourier component of the electric field, 

S_,(t) = <^(0) • E,(t)> . (1) 

k k k 


Hie notation means the following. E(x, t) = £* E+(t) exp(ik • x) is 

k k 

the electric field produced by N positively and N negatively charged 
rods of charge per unit length (± efl). The "particles" (rods) are 
aligned parallel to a uniform, constant magnetic field B = Be 

z 

and are located by giving their x and y coordinates. End effects 
of the rods are neglected. The velocities of the rods are given by 
the local cE x B/B drift. The electrostatic field E is self- 
ccnsistently determined through Poisson's equation. Periodic bound- 
ary conditions in the x and y direction are assumed, and the allowed 
values of k are k = 2n(n x , n^)/L, where n x , n y are any integers, not 
both zero. The symbol ( ) means an ensemble average over the t = 0 
positions of the particles, assuming they are distributed according 
to the canonical ensemble of Gibbs. 

Two 1 ' 7 theories of the temporal evolution of S_ > (t) lead, 

k 

making different dynamical assumptions, to strikingly different 



k 


predictions. The theory of Taylor and McNamara, 1 based on an assump- 
tion that the electric field seen by a particle obeys a jointly 
normal probability distribution at successive instants of time, leads 
to a behavior 


S™(t) = < Je | > exp[-k 2 R(t)] . (2) 

k k 

Here R(t) is a monotonically increasing function of t which obeys 
d 2 R(T)/dT~ = -9V(R)/3R, where V(R) is given by 

-(c 2 /4B 2 ) 2* <| E -»! 2 > (1 - exp(-2k 2 R))/k 2 , with R(0) = 0, dR(0)/dT = 0. 
k k v 

As t ■* », R(t) -» Dt/2, where D is, up to a factor of J£ f the Taylor- 
McNamara diffusion coefficient. 5 On the other hand, Taylor and 
Thompson 7 have, by means of a variant of the random phase approxima- 
tion, derived an expression 


s^(t) = ( |eJ 2 > cos(n^ '-) , (3) 

k k k 


an undamped sinusoidal oscillation for which 10 H. 2 -* 0(k ) for k small, 

k 

and is a more complicated function of k for k finite, but which is 
always real. The undamped oscillatory behavior of Eq. ( 3 ) contrasts 
sharply with the monotonically damped behavior of Eq. (2), and our 
purpose here is to determine via a numerical simulation of the system 
which, if either, prediction is correct. 
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II. SIMUIATION 


A program to follow the dynamical evolution of U000 rods 
advances the particle positions by the local I x B drift at each time 
step. E is determined through fbisson's equation by means of fast 
Fourier transforms and particle-in-cell (P.I.C.) methods, employing 
area weighting. 

Initial positions appropriate to thermal equilibrium are 
something of a problem since they are not capable of being produced 
by random number generators in any elementary way. We determined 
them from an unmagnetized two-dimensional simulation code which 
employed Newton’s laws of motion to advance the particle positions. 

This program was allowed to run, starting from random initial coordi- 
nates, through a few previously-determined 11 thermal relaxation times. 

The spectrum (e^ p advanced from the l/k^ form characteristic of ran- 

k / 

dom loading to the shape shown in Fig. 1. These final coordinates 

were then used as initial data for the guiding-center runs. 

Eight runs were carried out. Conservation of energy was main- 
tained to three per cent throughout . Autocorrelations were computed 
from three separate time origins, giving us an "ensemble" of twenty- 

four values of E_,(0) * E_^ ( v ) for each value of t. 

k k ^ 

Typical values of S ,(t) for different k modes are shown in , 

k 

Figs. (2a), (2b), and (2c). The crosses are the measured experimental 
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points, and the curves labeled "T & M" are Eq. (2), normalized to 
its initial value, t is measured in the appropriate dimensionless 
units. 12 

The curves in Figs. 2 are typical ensemble averages. A com- 
parison of eight individual runs with the ensemble average of twenty- 
four is made in Fig. 3. (The ensemble averages are crosses.) The 
large fluctuations about the ensemble average are characteristic of 
the runs and are believed to be genuine and not due to numerical 
error. Such accuracy checks as we have on the program (energy con- 
servation, independent tests of the Boisson solver, tests of the 
particle advancement algorithm in the analytically soluble case of 
two particles in a box) indicate substantially smaller accumulated 
errors by, say, t = 20 than the fluctuations indicated in Fig. 3. 

Theorems on the smallness of time-dependent fluctuation quantities 
about their ensemble averages are conspicuously absent in the theory 
of the guiding-center plasma. Persuasive theoretical arguments 

I 

exist to indicate that such theorems are not to be expected, unlike 
the more usual case of the plasma in which no guiding-center approxi- 
mation is made. 

In Figs, (4a) and (4b) a comparison is made between the ejqperi- 

ment&l values of S^(t) and Eqs. (2) and (3), normalized to their 
k 

initial values. Equation (2) is the curve labeled "T & m" and Eq. (3) 

( 

is labeled "T & T". It is clear that while the quantitative fit of j 

neither theory is perfect, the phenomena are correctly described 
qualitatively much better by Eq. (2) than by Eq. (3). 
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In Pigs. (5a) and (5b) the behavior of the autocorrelation is 

shown for various k- values, and a plot of k R (t) vs k L /4tt is 

given for various values of k at t = 20. Figure (5b) indicates a 

dashing of the autocorrelation which varies correctly with increasing 
2 

k , but indicates also a systematically weaker damping than the 
theory predicts. 
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III. DISCUSSION 


We have provided what is believed to be a rather sharp experi- 
mental differentiation between two competing theories of the temporal 
behavior of electric field autocorrelations in a two-dimensional 
guiding-center plasma. We have shown the qualitative correctness of 
the theory of Taylor and McNamara. 1 A systematically weaker decay 
than predicted by theory has been observed, a phenomenon for which 
at present no satisfactory theoretical explanation exists. Probably 
it is connected with the finiteness of the discreteness parameter 
(n Q \£)” , since the essential approximation in reference 1 seems to 
be the neglect of correlations between single-particle orbits. 
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FIGURE CAPTIONS 


Electric field fluctuation spectrum vs wavenumber, normalized 
to initial value. Solid line is theoretical thermal equili- 
brium curve derived from Debye-Huekel theory. Crosses are 
time-averaged values for an unmagnetized simulation plasma 
in two dimensions, where average is carried out over several 
relaxation times. Circles are single-time values and are 
present to indicate the level of typical fluctuations. 

Autocorrelations S_,(t) vs t for three different modes, 
k 

Crosses are measured experimental points and curve labeled 

"EXP" is a fit to these points. Curve labeled "T & M" is 

Eq. ( 2 ). Modes are labeled by values of (n , n ). 

x y 

(2a): (n x , n y ) » (4, 4); (2b): (n x , n y ) » (0, 6); 

(2c):(n x , n y ) = (6, 6). 

Autocorrelations S.(r) vs r. Eight individual runs are 
k 

compared with ensemble average values over a sample of 

twenty- four. Crosses are the ensemble average . 

Comparison of the results of experiment (crosses) for S.(t) 

K 

with theoretical values given by Eq. (2) (labeled "T & M") 

and Eq. (3) (labeled "T &T"). (4a) is for (n , n } = (4, 4) 

x y 

and (4b) is for (n . n } » (0, 6). 

* jr 
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Fig. 5 


(5a) is S.(t) vs t for vario»?K k values. (5b) is 

5 ? 

P / ^ p 

k R (t = 20) vs n + n . The solid curve is the theoreti- 
x y 

cal prediction of Taylor mid McNamara. 
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Fig. 2b 
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Fig. 2c 
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Pig. 4b 









